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Abstract 



We study the dynamics of classical and quantum systems undergoing a con- 
tinuous measurement of position by schematizing the measurement apparatus 
with an infinite set of harmonic oscillators at finite temperature linearly cou- 
pled to the measured system. Selective and nonselective measurement pro- 
cesses are then introduced according to a selection of or an average over all pos- 
sible initial configurations of the measurement apparatus. At quantum level, 
the selective processes are described by a nonlinear stochastic Schrodinger 
equation whose solutions evolve into properly defined coherent states in the 
case of linear systems. For arbitrary measured systems, classical behavior is 
always recovered in the macroscopic limit. 
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I. INTRODUCTION 



A fundamental problem in quantum mechanics is the relationship between the states in 
the Hilbert space of a quantum system and the states in the phase space of the corresponding 
classical system. This is particularly evident in the case of a superposition of states which are 
individually mapped in the macroscopic limit, formally ^ ^ 0, into distinguishable classical 
states 0,1 . 

One important step toward the solution of this problem has been made by recognizing 
that a system is never completely isolated by the external world. It has been argued that 
an external environment can, after a transient whose duration presumably depends on the 
coupling strength, drive the totality of the admissible states of the Hilbert space into those 
having macroscopic limit 

Among all the conceivable situations which require the interaction with an environment, 
a peculiar role is played by the measurement processes. Indeed, whenever any physical 
property of a system is investigated, an unavoidable coupling with the degrees of freedom of 
the measurement apparatus must be invoked. Taking into account these external degrees of 
freedom naturally provides a generalization of the von Neumann postulate to continuous 
measurements f^-^. 

There exists a basic difference between a general environment and one schematizing a 
measurement process. According to the Copenhagen interpretation, classical behavior of 
the measurement apparatus has to be assumed before the information is registered by the 
observer |]5|, p!0|JTT[| . This requirement, although establishing a classical connection between 



observer and meter, does not imply a classical behavior of the observed system since its 
interaction with the meter fully preserves quantum features. This aspect may be emphasized 
by comparing the coupling of the same measurement apparatus, classically controlled by the 
observer, to classical or quantum systems. 

This paper is devoted to establishing the relationship between the dynamics of classical 
and quantum systems under the influence of a measurement process. We model the effect of 
the measurement by allowing the measured system to linearly interact with an infinite set 
of harmonic oscillators. The interaction occurs in the configuration space thus representing 
a direct measurement of position. Our approach shouldn't be understood as the modelling 
of a particular device, e.g., a cloud chamber in which the oscillators represent molecules 
distributed in a medium |0. Our oscillators are the modes of an interaction field which 



reproduces the main features expected in a measurement process, e.g., the wavefunction 



collapse, without ad hoc rules. Differently from other abstract models |]T3[, our approach 
has, however, a certain degree of realism since allows us to operatively define the reading of 
the measurement results. 

The oscillators representing the measurement device are chosen at thermal equilibrium 
with temperature T and continuously distributed in frequency with a proper density. These 
conditions ensure that, for a chosen initial configuration of the oscillators, the classical 
measured system is described by a Markovian Langevin equation with white noise and a 
certain relaxation time 7~^. The constant 7, which fixes the magnitude of the density of 
oscillators, represents the strength of the measurement process. In the quantum case, a 
further characteristic time of the measurement apparatus, h/ksT, arises. The requirement 
of classical behavior of the meter with respect to the observer implies a high temperature 
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condition to hold so that this thermal fluctuation time must be much shorter than the 
relaxation one. As a consequence, the measured quantum system is described by a stochastic 
nonlinear Schrodinger equation p^-p . 



Details of the derivation of the stochastic equations describing classical and quantum 
systems during a measurement process are given in Sections II and III, respectively. We call 
these processes selective since they correspond to a single measurement act with initial con- 
ditions of the meter selected among those compatible with the assumed thermal equilibrium. 
Alternatively, one can also consider nonselective measurement processes corresponding to an 
average over all initial configurations of the meter. In this case, the dynamics of classical and 
quantum measured systems is described by a Fokker-Plank equation and a trace-preserving 
positive master equation, respectively. The experimental results obtained by repeating a 
measurement on a system always in the same initial state or performing the same measure- 
ment on an ensemble of equally prepared independent systems is directly comparable with 
the solution of these nonselective equations. Equivalent results are predicted by averaging 
the solutions of the selective equations, classical or quantum, over the corresponding noise 
realizations. 

In Section IV we show how to infer the measurement results by the reading of an appro- 
priate pointer. Since the oscillators representing the meter are classical with respect to the 
observer, the pointer can be defined in terms of the coordinates of these oscillators which, 
in turn, reflect the status of the measured system. 

As for closed quantum systems, also in presence of a measurement process there exists 
a class of states, namely the coherent ones, which admits, in a proper sense, the ^ — 
limit. These coherent states, explicitly built in Section V, are Gaussian localized states in 



the co-moving frame of a measured linear system [|1^,|T^,0. In Section VI we show that for 
linear systems the solutions of the stochastic Schrodinger equation converge to a coherent 
state localized around a point in the phase space which moves according to a Langevin- 
like equation. This equation reduces to the classical one for h ^ 0. The convergence 
into this coherent state occurs in a timescale the geometric mean of the two 

characteristic times associated to the measurement apparatus. This time diverges for an 
unmeasured system (7 — > 0) and vanishes in the macroscopic limit. 

In the case of nonlinear systems, the phase-space localization through convergence into 
a coherent state becomes the leading dynamical process as ^ — > 0. This is sufficient to 
demonstrate that classical behavior is always recovered in the macroscopic limit, avoiding 
any paradoxical quantum feature. 



II. CLASSICAL SYSTEMS 



Let us consider a system described by the Hamiltonian 

2 

H{p,q,t) = ^ + V{q,t) (2.1) 

and suppose that we want to measure its position q {q & R, for simplicity). We schematize 
the measurement process by the interaction of the system ( p^.l| ) with a set of particles of 
mass M and canonical coordinates {Pn, Qn) via a harmonic potential which, for evident 
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physical reasons, must depend on the relative distances Q„ — q. The Hamiltonian for the 
total system is then taken as 

Htot = Hip,q,t) + HmiP,Q-q), (2.2) 

where 



HUP,Q-q)=J2 



2M 



2.3) 



represents the measurement apparatus linearly coupled to the measured system. Here and 
in the following, Q and P are a shortening for the whole set {Qn} and {Pn}- Note that the 
Hamiltonian (|2.3|) can be interpreted as that of a set of harmonic oscillators with equilibrium 
positions Qn = q- Our model is thus reduced to the exactly solvable problem of a system 
interacting with a bath of harmonic oscillators |2^ but with the difference that the inter- 



action potential is invariant under space translations. The importance of this invariance in 
avoiding the appearance of infinite renormalization potentials has been recently underlined 
in |21] in the framework of the classical/quantum Brownian motion. We stress that the long 
range nature of the quadratic potential in ( |2.3|) prevents us from considering a situation in 
which the measured system does not interact with the meter so that our model cannot take 
into account the switching-on of the measurement. 

At classical level, the Newton equation for each harmonic oscillator can be explicitly 
solved in terms of the function q(t) and the values of the coordinates Q'„ = Qn(t'), P^ = 
Pn{t'), and q' = q{t') at the initial time t' 

Quit) = {Q'n - q') cosK(t - t')] + sin[cu„(t - t')] 

+q{t) - [ ds cos[uJn{t - s)] q{s) . (2.4) 
Jt' 

We will also assume p{t') = p'. When this solution is inserted in the Newton equation for 
the measured system, the following equation for q(t) is obtained 

mq{t) + r ds r{t - s) q{s) + dqV{q{t),t) = U{t), (2.5) 
Jt' 



where 



and 



T{t-s)=J2Mujlcos[uJnit-s)] (2.6) 



(2.7) 



In the limit of an infinite number of oscillators, if the corresponding initial conditions Q' 
and P' are a realization of a stochastic process in the phase space, H(t) is a realization of a 
stochastic process in time. If, as we will suppose, the oscillators are at thermal equilibrium 
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with temperature T, the initial conditions to be considered are typical realizations of the 



stochastic process corresponding to the equilibrium Gibbs measure In this case, the 

following statistical properties for 11 (t) hold 



n(t) = 0, U{t)U{s) =kBTT{t-s), (2.8) 

where 

_ / dP'dQ' ■ ■ ■ e-^™(^''«'-^')/'=s^ 



JdP'dQ' e-Hm(P',Q'-g')/kBT 



(2.9) 



Note that the initial conditions corresponding to the definition (|2.9|) respect the space trans- 
lation symmetry of the Hamiltonian Hm thus implying, in agreement with the long range 
nature of the quadratic potential, a correlation between the coordinates Q' and q'. 

The friction term in the stochastic differential equation (|2.5| ) contains memory effects 
which are an unessential complication in our context. A Markovian evolution can be obtained 
by choosing an appropriate continuous distribution of the frequencies {c<j„}- For a frequency 
density 

du ttMuj^ V v ; 

where 6{x) is 1 for x > and otherwise, we get 

r°° , dN r, ^ , sin[f2(t — s)l ^, , ^ 

T{t-s) = j^ du—Muj^ cos[a;(t - s)] = 2m-f - 2m-^6{t - s). (2.11) 

The approximation holds for t — s^ Q^^. For ^ r, where r is the fastest time scale at 
which the measured system evolves, 11 (t) can be approximated with a white noise and Eq. 
TSl) rewritten as 



dp{t) = - [7p(t) + dgV{q{t),t)] dt + ^2m-fkBTdw{t) (2.12) 

dq{t) = ^dt. (2.13) 

m 

Here, we introduced the Wiener process dw{t ) = {2m-fkBT)-^/^U{t)dt having zero average, 
dw{t) = 0, and standard scaling, dw{t)dw{t) = dt. Note that the friction coefficient 7 and 
the temperature T completely define the fluctuation and dissipation phenomena induced by 
the interaction with the measurement apparatus. 

Equations (|2.12|) and (|2.13|) describe the evolution of the system during a selective mea- 
surement, i.e., a measurement in which a realization of the stochastic process p{t), q{t) is 
selected according to certain initial conditions of the measurement apparatus. Alternatively, 
one can consider a nonselective measurement corresponding to an average over all possible 
realizations of the stochastic process. In this case, the measured system is described by a 
probability density W{p,q,t) which is determined by the Fokker-Plank equation ||23| asso- 
ciated to (1^ ) and (|]T§ 



dtW{p,q,t) 



-—dq + dqV{q, t)dp + dp (7^ + m-fkBTdp) 
m 



W{p,q,t) (2.14) 
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with initial conditions W{p,q,t') = 5{p — p')S{q — q'). The probabihty density W{p,q,t) 
allows us to directly evaluate averages of any function of p{t) and q{t). In particular, for the 
position we have the average value 



dpdq W{p, q, t) q 



(2.15) 



with a variance 



Aq\t) 



dpdqW{p,q,t) q — q{t) = q{ty — q{t) 



(2.16) 



For 7 — i> the effect of the measurement vanishes. In this case, Eq. (|2.14|) becomes the Liou- 
ville equation for the isolated system ( |2.1| ), the average value ( p^.l5D gives the corresponding 
time dependent solution, and the variance ( |2.16|) vanishes. 



III. QUANTUM SYSTEMS 



At quantum level, the measured system is conveniently described by a reduced density 
matrix obtained by tracing out the coordinates of the measurement apparatus in the total 
density matrix 



^'(^1, 92, t) = J dQ QtoMu Q, 12, Q, t). 



(3.1) 



We assume that at the initial time t' the oscillators of the measurement apparatus are at 
thermal equilibrium with temperature T and the total density matrix is factorized as 



QtoM'i, Q'l, Q2, t') = g{q[, q2, t')gmiSQ[, SQ'^, t' 



(3.2) 



where 



MuJr. 



■ tanh 



\\ -kTi \2kBT 



X exp 



/2 

2n 



2n 



2h \taj\\i{huLin/kBT) sm\i{huJn/kBT) 



(3.3) 



with 5Q'i = Q'l — Qeq and 6Q2 = Q2 — Q'eq- The requirement of translational and time 
reversal invariance implies pl| 



Q'e 



<i'i + q'2 



eq 



(3.4) 



Analogously to the classical case, this choice corresponds to an initial condition ( p.2|) fac- 
torized but correlated. 

At a later time t, the reduced density matrix is obtained through a Green function 



Qi'li,'l2,t) = J dq[dq'2 G{qi, q2,t; q[, q2,t') giq[,q2,t') 



(3.5) 



whose path-integral representation 

Giquq2,t;q[,q2,t') = J rf[pi]rf[gi]^|;*, J d[p2]d[q2]ll'l, 

X exp|^5'[pi,gi] - ^S[p2,q2] - Z[pi, qi,p2, q2]j (3.6) 
is the free evolution of the measured system modified by the influence functional 
exp{-Z[pi,gi,p2,g2]} = j dQ j dQ[dQ'^J d[Pi]d[Qi]^l^, J d[P2]d[Q2]'^l,, 

X exp!^^Sm[Pi,Qi-qi] - ^S^[P2,Q2-q2]]Qm{SQ'i,SQ'„t'). (3.7) 

Here, S[p,q] and Sm[P,Q — q] are the classical actions corresponding to the Hamiltonians 
H and Hm, respectively 

S\p,q]= [ ds\p{s)q{s) - H(j),q,s)] (3.8) 
Jt' 

S4P,Q-q]= f ds[P{s)Q{s)-HUP,Q-q,s)]. (3.9) 
Jt' 

The functional measure with boundary conditions q{t') = q' and q{t) = q is obtained by 
slicing the interval [t', t] at times t^'^^ = t' + {t — t')n/N, n = 1, . . . , N, and taking the N ^ oo 
limit. All time integrals can be approximated with sums 

ldsf{p,q,s) = Y: I dsf{p,q,s)c^-—Y.f{p^-\q^-\t^-^), (3.10) 

where p*^"^ = p{t^"^) and g*^"^ = q{t^^^), and 

N ^„(n) N-1 

i'^"-ii'Lni^n/^'"'- (3.11) 



n=l n=l 



Analogous relations hold for d[P]d[Q]Q', . 

The influence functional (|3.7|) contains only Gaussian integrals and can be evaluated 
exactly. The result is [E^ 



Z[pi,qi,P2,q2] = f ds f du [qi{s) - q2{s)] 
Jt' Jt' 



where 



X ^ A(s - u) [q,{u) - q2{u)] + ^r(s - u) \p^{u) + P2{u)] \ , (3.12) 



Hs-u) = ^2^^„coth I I cos[uJn{s-u)\ (3.13) 
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and r(s— m), defined by Eq. ( |2.6| ), are named fluctuation and dissipation kernels, respectively 
20| . The double time integral which appears in Eq. ( 3.12|) is responsible for memory effects 



which break the semi-group property of the evolution of the measured system, i.e., 

^(gi, g2, t; q[, t') = J dq'ldq'^ G{qi, ga, t; g", q'^, t") G{q'l, q'^, t"; q[, q'^, t') (3.14) 

for t' < t" < t. Two are the sources of this non-Markovian behavior. In the dissipation kernel 
r, the origin of memory effects is classical and Markovian behavior is obtained if one assumes 
the frequency distribution ( p^.lOD with -C r. In the fluctuation kernel A, the assumption 
of these conditions does not remove the memory effects due to the quantum behavior of 
the oscillators in the thermal bath. However, in the framework of a theory of measurement 
processes, the measurement apparatus, which is an interface between the observer (classical) 
and the measured system (quantum), must have classical behavior with respect to the former 
in order to avoid paradoxical features |0. In the present model, the nature of the interaction 
between the measurement apparatus and the external environment is fixed by the ratios 
hun/kBT SO that the requirement of classical behavior of the measurement apparatus with 
respect to the observer imposes the condition fc^T ^ hQ. 

For h/ksT -C ^ r and frequency distribution ( p.lO| ), Eq. ( |3.12| ) can be approxi- 
mated (see Appendix A) with 



Z\pi,qi,P2,q2] = J^^ rfsj 



j m-fknT ^^^^^^ _ ^^^^^^2 ^ ^ ^^^^^^ _ ^^^^-^^ ^^^^-^ +P2is) 



(3.15) 



A differential equation for the reduced density matrix operator can be then derived with 
standard methods. According to Eq. ( p.l5|) and using the path-integral representation of 
G, we have 

QiQi, 12, t + dt) = J dq[dq2 G{qi, q2, t + dt; q[, q'^, t) g{q[, q^, t) 

dpi , f dp2 , (i , , . i , , . 



+ 



'2h 



ill - (12) {Pl + P2) 



dt \Q{q'i,q'2,t) ■ 



(3.16) 



By using the identity {q\p) = {2ttK) exp(ipq/fi) and expanding the exponential containing 
the infinitesimal time dt, the above equation can be cast in the form 



{qi\git + dt)\q2) = {qi\git)\q2) + {qi\-^g{t)\q2) , 



where 



^^g{t) = -^[H{p,q,t),g{t) 



-^2 — % ^(^)]] - ^ ^(^)}] 



(3.17) 



(3.18) 



Note that d[TTg{t)\/ dt = so that we can assume Tr^(t) = 1. 
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Equation (|3.18|) describes the evolution of the measured system with initial conditions of 
the oscillator system statistically distributed according to Eq. ( p.3|) . This is a nonselective 
measurement process to be compared with the classical one described by Eq. ( |2.14| ). 

The quantum-classical correspondence of nonselective measurements can be extended to 
selective processes. Introducing 



A{p{t),q{t))- 
B{p{t)) - 

Eq. (|3.15|) can be rewritten as 



'2m'ykBT 



7 



SmkBT 



Pit) 



7 



SmkBT 



Pit), 



(3.19) 
(3.20) 



ds 



ds 



«7 
^7 



1 



iqiPi - q2P2) + -AiAl + -A2A; - A^A; - -(5i - B: 



2 ' 2 

1 ... 1 ... 1, . ,2 1 



2^(giPi - q2P2) + -AiAl + -A2A; - -{Ai - a)' - -{A; - a*)' + aa* 



-A^a* - Ala + ^(Ai -a-Al + a*)^ - ^(fii ~h-B2 + hf 



(3.21) 



where Ai and Bi stand for A{pi,qi) and B{pi), i = 1,2, and a{t) and b{t) are two arbitrary 
functions, respectively complex and real. The coupling between the components 1 and 2 of 
the system coordinates given by the last two squares of the exponential can be eliminated 
in terms of two functional integrations over white real noises by means of the identities 



exp 



1 /■* 

— / ds {Ai-a- AI + a*)^ 

2 «/ 



d[^] exp 



ds 



(Ai - a)' + {Al - ay - {A, - a)^ - {Al - a*)^ 



(3.22) 



and 



exp 



I f ds {B,-b-B2 + hf 



d[rj\ exp 



ds 



(5i - 6)2 + {B2 - h)^ + i{B^ - b)ri + i{B2 - b)ri 



(3.23) 



Note that the above functional integration measures are Gaussian 



din = lim n d^ 

n=l 



It-t' 

2txN 



exp 



it-t')i ^^y 

2N 



so that 



m = J d[^] m = 0, ams) = J d[^] m a^) = sit - s) 



(3.24) 



(3.25) 



and analogously for 77. The two-particle Green function ( p.6| ) can be then rewritten in terms 
of a couple of one-particle Green functions 
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GiQi,Q2,t;q[,q2,t') = J d[^]d[ri] G^^^^{qi,t; q[,t') G^^^^{q2,t-q'^,t' 



(3.26) 



where 



d[v\dW\f,t' exp 2h It' 



+ / ds 



^AA* -^{A-aY - ^aa* + Aa* + (A - a)^ + {B - hf ± i{B - 6)r/] |. (3.27) 

Assuming that the system is initially in a pure state, i.e., git') = \4'(t')){ip(t')\, at a later 
time t the reduced density matrix operator is expressed as a functional integral over pure 
states 



(3.28) 



obtained propagating \i'{t')) with Gj^^j 

{<l\i^Mit))= fdq' G|,](g,t;g',t') (qW))- 



(3.29) 



An evolution equation for the states |V'[^,,](^)) can be obtained by writing the exphcit 
form of the propagators G^^^] between the times t and t + dt 



{Qli^^it + dt)) 



i G%Aq,t + dt-q\t) {q'm)) 



-]^A{p, q)A{p, q)* - ^[Aip, q) - a{t)f - \a{t)a{ty + A{p, q)a{ty 



+ [A(p, q) - a{t)m + [B{p) - h{t)f ± i[B{p) - h{t)Ut) 



dt 



x(gVL(^))- 



(3.30) 



By using the identity {q\p) = (27r/i)~^/^ exp{ipq/h) and expanding the exponential containing 
the Wiener processes dw^{t) = ^{t)dt and dwrjit) = rj{t)dt according to the Ito rule [^, we 
get the following stochastic differential equation [^] 



7 

H{p,q,t) + -{pq + 



2 
1 

+ 2 



+ a{tya{t) - 2a{tyA\ |V^|^](t)) dt+\A- a{t)] |^f^„i(t)) dw^{t) 



B - h{t)\ l^l^j(t)) dt±i[B- h{t)\ |^±^,(t)) dw,{t). 
The normalization condition for the reduced density matrix operator 

d[C\dVl] =1 



T,Q{t)= dmr,] (^[^,](t)|^J,](t)) 



(3.31) 



(3.32) 
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is satisfied by imposing {4'[^ri]i^)\4'^r]]i^)) = 1- This fixes the arbitrary functions a(t) and 
b(t). Indeed, the requirement that the Ito differential 



- a{t) + it _ a*(t)|^J^](t)) dw^it) 
+2^(^[-^j(t)|5-6(t)|^+^](t)) dw.it) 



(3.33) 



vanishes, imphes 



l2m'ykBT 



7 



7 



Y SmkBT 



(3.34) 
(3.35) 



The appearance of two different Green functions G^^^ in Eq. 



).28| ) possibly introduces 

a violation of positivity in g{t). This unphysical property is reflected in the anomalous 
definition of expectation values of Hermitian operators, e.g., ('^[^^](^)|^|'?/'[^^](^)), which may 
be complex. The problem is mathematically related to the presence of the B terms in Eq. 
( p.21| ). However, due to the high temperature condition h/ksT -C r we have B <^ A (see 
Appendix A for details) and the last square in Eq. ( |3.21|) can be neglected with respect to 
the last but one. In this case, Eq. (|3.18| ) becomes 



rwykBT 



[q, [Q,m]]-7;^% {p,m}] 



7 



2h 



IQmkBT 



+ 



Q{t) 



1 

+ 2 



(3.36) 



This equation is of Lindblad class and provides a (completely) positive evolution of Q{t) 
The reduced density matrix operator can be decomposed in terms of a single state 
associated to the Green function G[g] obtained by neglecting the B terms in |'0fc„i(^)) and 



Gj^^], respectively. Note the disappearance of the rj noise. Equation ( p.31|) becomes the 
norm-preserving stochastic Schrodinger equation 



d\i^m{t)) 



I 
n 

1 r 
2 



7 



+ 



](t)) dt 



A^A + a(t)*a(t) - 2a(t)*ij \i'[^]{t)) dt + lA - a{t)] \^[^]{t)) dw^{t) (3.37) 



with a(t) = (^[^](t)|i|VK](t)). 

Equation (|3.37|) describes the evolution of the measured system for a realization of the 
stochastic processes \ip[^]{t)). This is a selective measurement process which is related the 
nonselective one by the relationship 



Pit) = l%](t))(%](t)l = / d[^] |%](t))(%](t)| 



(3.38) 
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The quantum expectation values of the observables, e.g., q(t) = {■ip[Q{t)\q\ip[Q(t)) , are 
stochastic processes with average value 



and variance 



Aq'it) = jd[i] (%](t)|[g-g(t)]'|%](t)) 



(3.39) 



(3.40) 



According to Eq. (|3.38| ), these average quantities can be also directly computed by consid- 
ering the nonselective measurement process described by Eq. ( p.36|) 



q(t) = Tr 



dpdq W{p, q, t) q 



(3.41) 



Aq\t) = Ti\m q-q{t) 



dpdq W{p, q,t) q — q{t) 



(3.42) 



The last expressions, formally identical to the classical Eqs. ( |2.15| ) and ( |2.16| ), have been 
obtained by introducing the Wigner function W{p,q,t) through the relation g{q,q,t) = 
J dp W{p,q,t). For 7 — >• 0, the effect of the measurement vanishes. In this case, Eq. ( ^.37[ ) 
becomes the Schrodinger equation for the isolated system and the variance ( p.42|) reduces 
to the standard quantum mechanical expression. 

Finally, we show that the von Neumann collapse theory is recovered for 'j 00 and 
T 00 with T7~^ constant. In this limit, we must identify the shortest time scale r 
of the classical system with 7"^ and the condition T7~^ constant allows the inequality 
h/kBT -C fl~^ <^ r to be always satisfied. At time t = t' + r, Eqs. (|3.36|) and ( p.37|) give 



g(qi,q2,t' + r) ~ exp 
ijj^/ (g, t' + r) ~ exp 



--Kr(gi - ga)^ 



'KT \ q — q 



g{quq2,t' 
2 



2v^, 



exp f ^ I ^iq,t') 



7 



(3.43) 
(3.44) 
0, Eqs. 



where n = 2m'ykBT/h'^ , ^' = ^{t'), and q' = q{t'). For T'j~^ constant and r 
( p.43| ) and (|3.44|) provide an instantaneous diagonalization of the reduced density matrix and 
an instantaneous collapse of the wavefunction into the eigenfunction of q corresponding to 
the eigenvalue jl^fn.. The normalization factor in Eq. ( p.44|) is such that the quantum 
expectation values at time t' + r, when averaged over the noise realizations coincide with 
the quantum expectation values at time t' . For instance, we have 



(V^e(t' + r)|g|^5,(t' + r)) 



dq q |^/'g/(g,t' + r)|^ 



'M2 



dq q \ij{q,t')\ 



(3.45) 



which is the result expected on the basis of the von Neumann postulate: in a selective 
measurement of position at time t', the probability that the state {ipit')) collapses into the 
eigenstate \q) is \{q\ip{t'))\'^. 



12 



IV. MEASUREMENT RESULTS 



In the previous Sections, we have seen how the evolution of a system, classical or quan- 
tum, is influenced by coupling its coordinate to those of infinitely many linear oscillators. 
We called this process a measurement of position in agreement with the fact that in a proper 
limit the von Neumann collapse theory can be recovered from the resulting equations. Now, 
we specify how the properties of the measured system can be operatively read by the ob- 
server. 

From the point of view of the observer the linear oscillators of the meter always have 
classical features so that their coordinates can be directly taken as pointers of the meter 
itself. Consider, for example, a pointer whose value R{t) is defined as 

R{t) = I dsY.Ut~s)Qn{s). (4.1) 

During a single measurement, i.e., in a selective process, R{t) is a stochastic variable. We 
would like to choose the response functions C^nit — s) in order that the statistical properties 
of R{t) over the ensemble of all possible measurements coincide or, at least, allow us to 
recover the nonselective properties of the measured system. For instance, we could ask that 
the average measurement result 

R(t)= fdsJ2Ut- s)Q4^ (4.2) 

and its variance 

AR\t) = JdsJ2 Ut - s) [q:;(^ - Q:;(s)'] , (4.3) 



correspond to the quantities ( ^.15| ) and ( |2.1(j| ) in the case of a classical system or those ( |3.41D 
and ( |3.42| ) in the case a quantum one. 



In the classical case, Eq. ( |2.4| ) provides an explicit expression of the oscillator coordinates. 
By using the property <^ r, we see that for ~ f2 we can find a period much 
shorter than the fastest classical time and much longer than the inverse of the oscillator 
frequency so that the average of Q„.(t) in this period coincides with q{t). With the choice 
Cn(t) = Aexp(— At) 5(^„-f7 and using the definition ( p.9|) we then have 



Rit) = q(t) (4.4) 

and 

Ai?'(t) = Ag2(t)+£2_ (45) 

The pointer variance is the sum of the variance Ag^(t) of the measured system and the 
resolution of the measurement apparatus 

2 kBr_ 
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The term i"^ represents a systematic error of the measurement and can, in principle, be 
subtracted. 

The above results can be derived in an alternative way. Consider the general definition 
of the moments 



Qn{t) = / dpdqdPdQ Wtot{p,q,P,Q,t) Qr 



(4.7) 



Qn{ty 



dpdqdPdQ Wtot{P: q, P, Q, t) Ql 



(4.8) 



Here, Wtot{p,q, P,Q,t) is the probability density solution of the Liouville equation for the 
total system with initial conditions 



WUp,q,P,Q,t') = 6{p - p')6{q - q') 



-Hm{P',Q'-q')/kBT 



JdQ'dP' e-^-(-f"'«'-9')ABT- 



(4.9) 



The oscillators with frequency cu^ ~ i7 approach the thermal equilibrium around the instan- 
taneous value of the measured coordinate on a time scale much shorter than a characteristic 
period with ^ ^ r. For these oscillators the time average of the moments 
( [4.7[ ) and ( [4.8| ) over a period A^^ can be approximated by inserting the following adiabatic 



expression for the total probability density in the same Eqs (^4.7|) and (^]8| 

^-H,n(P,Q-qit))/kBT 

Wtot{p,q,P,Q,t)^W{p,q,t) 



(4.10) 



JdQdP e--f^-(-P.<3-'?W)As^' 
where W{p,q,t) is the solution of Eq. ( |2.14|) . Equations ([4.4]) and ( |4.5| ) are then obtained 



by evaluating the Gaussian integrals in (|4.7| ) and (^ 



The last approach formally applies also at quantum level by interpreting the Ws as 
Wigner functions. The adiabatic approximation ( [4.10 ) becomes 



Wtotip, q, P, Q, t) ~ W{p, q, t) WUP, Q - q{t)), 
where W{p,q,t) is the Wigner function associated to g(t) and 

tanh (^^) ( / huj \ I 

WUP, Q - qit)) = n exp \ - tanh ' ^""^ ^ 



(4.11) 



2kBT 



Moo. 



hMuj„ 



h 



- {Qn - q{t)f 



(4.12) 



is the Wigner function associated to the density matrix ( ^l3|) with Q^q = q(t). Here, q(t) is 
the quantum expectation value of q in the state |?/'[g](t)). Due to the condition ksT ^ hfl, 
the above expression reduces to the classical distribution 



WUP,Q-qit)) = l[ 



UJr, 



2TxkBT 



exp 



(4.13) 



so that Eqs. 
apparatus. 



and 



still hold with the same resolution I of the measurement 
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V. THE ;i ^ LIMIT: DEFINITION OF COHERENT STATES 



The dynamics of a closed quantum system reduces to the classical dynamics in the h ^ 
limit only if the system is prepared in appropriate states. The coherent states defined as 
the ground state of the displaced harmonic oscillator 



bV) 



(5.1) 



where |0o) is the ground state of the undisplaced oscillator, are a well known example p8| . 
These states provide a convenient representation for studying the h —>■ limit regardless of 
the nature of the Hamiltonian which may not preserve their form |^ . 

In the case of the measurement model discussed here, it is possible to find states which 
are localized and stationary in the comoving frame of the measured system and play the role 
of the ground state |0o) in Eq- ( |5.1| ). A first example of these states was given in [|1^,|15] for 
a free particle evolving according to the dissipationless Eq. ( |A6| ). A generalization valid in 



the case of a harmonic oscillator described by Eq. ( ^.37| ) has been recently provided in 
Here, we derive the expression of the coherent states for a general linear system with con- 
stant proper frequency undergoing the selective measurement process of Eq. ( ^.371 ). Then, 
we discuss the recovering of the classical limit in selective and nonselective measurement 
processes on arbitrary systems which are prepared in such states. 

During a selective measurement, the quantum system is described by a state \ip[^]{t)) 
which evolves in a specified rest-frame according to Eq. (|3.37 ). In analogy with Eq. (^.11 
we look for solutions of Eq. ( p. 37] ) of the form 



(5.2) 



where p{t) = {ip[^-]{t)\p\ip[^-]{t)) and q(t) = (■?/'[^](t)|g|'0[^](t)). The co-moving state |0) is 
assumed constant so that the solutions (|5.2|) depend on the noise C,(t) only through the 
expectation values p{t) and q{t) and the action (f{t). By inverting the transformation 
and imposing that the change of |0) in a time dt vanishes, we get 



,i[vp{t)+rf^(t)] 



^j.[p{t)+dp{t)]q ^l.[q(t)+dq{t)]p 



|%](t)) + rf|^[^](t)) = 



(5.3) 



The differential d\ip\^^{t)) is given by Eq. ( p.37|) . The same Eq. ( p.37|) allows the evaluation 
of 



dp{t) = - hp{t) + (^[5](t)|9,V^(g,t)|V^[g](t))l dt + 2^al^dw^{t) 



and 



dq{t) 



Pit) 



m 



dt + ( 2v^(tJ - ) dw^{t), 




where 



(5.4) 



(5.5) 



(5.6) 



and 
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0", 



pq 



77(%]WlM + gp|V^K]W) - (%](t)|p|V^K](t))(V^K](t)|g|%](t)) = ^{m + qm (5.7) 



are the constant variances associated to the states ( |5.2|) and k = 2m'^kBT /h^ . The expec- 
tation value of the force operator which appears in Eq. ( |5.4| ) can be expressed in terms of 
the co-moving state by a translation q ^ q + qif) 



(V^[5](t)|9,F(g,t)|^K](t)) = {<P\d,V{q + q{t),tm. 



(5.8) 



Finally, we write the differential of the action ip{t) in terms of two coefficients /i(t) and v{t) 
to be determined later 

dipit) = i2{t)dt + p{t)dw^{t). (5.9) 
By expanding the exponentials and using the Ito rule, Eq. (|5.3|) can be rewritten as 



i + Fdw^{t) + Gdt 



0, 



(5.10) 



which is equivalent to F\(j)) = and = 0. In general, the operators F and G will 

depend on time through the expectation values p{t) and q{t) and the action (f{t) so that 
these equations cannot be satisfied with a constant However, we can try to make F and 
G time independent with a proper choice of the coefficients fi{t) and z^(t). In the case of F, 
we have 



2i 



2i 



1 - —(Tpg ] q + —a J 



+ 



1 



2^, 



(5.11) 



and this becomes time independent with the choice 



it) = -p{t) (^2v^a, 
The corresponding equation F 

0(g) = {q\<P) = (2 

For z/(t) given by Eq. ( p.l2|) , the operator G is 

+ 



2 7 

" 2^^ 



(5.12) 



has the unique normalized solution 

-1/4 



1 - 



exp 



4a2 



(5.13) 



2m 
7 
2 



+ V^(g + q{t),t) - (0|a,F(g + g(t), t) |0)g - /i(t) 

1 



- 2Ka^galg 



(5.14) 



This can be made time independent with a proper choice of fi{t) only for linear systems. 
Assuming V{q,t) = Vo{t) + Vi(t)q + ]^mujQ(f with uiq constant, we have 



2-2 



V{q + q{t),t) - {<i)\d,V{q + g(t), t) |0)g = V{q{t),t) + -mu'^q 



(5.15) 
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so that by choosing 



the equation G\(j)) = becomes 



Pit) 



1 



2^ +V{q{t),t) + ^p{t)q{t), 



(5.16) 



2 ^2 



(5.17) 



In the position representation and using 0(g) given by ( |5.13|) , Eq. (|5.17|) is equivalent to the 
following two complex equations 



2m 



+ 



muo, 



2i „2 







i-2 1 2j^2 

2m 2al 



2 + 2 



1 - f 



^ - z^K = 



0"„ 



o 2 2 I -i- 2 



(5.18) 
(5.19) 



Equation ( ^.18|) gives two conditions for the determination of and 0"^^. The constant e 
can be then evaluated from the real part of Eq. ( |5.19| ) , the imaginary part being an identity. 
The solutions are 



7^ - uj^ + - ooi)^ + {2hK/my 



a. 



pq 



and 



(5.20) 
(5.21) 

(5.22) 



The variances and are always real and positive except for ksT/huJo <^ 1 which is, 



however, outside the range of validity of Eq. ( 3.37 ). According to the choices ( 5.12 ) and 
( ^.16| ), we finally have 



dLp{t) 



e + ^ + V{q{t),t) + jp{t)q{t) 



7 



dt-p{t)dq{t). 



(5.23) 



This allows for an interpretation of e in terms of a zero-point energy which adds to the 
classical renormalized Hamiltonian H{p{t),q{t),t) + ^p{t)q{t). 

In the 7 — i> hmit, we have cr^ = h/2mujQ, a'i = 0, and e = hujQ/2. The stationary state 



of Eq. ( ^.13|) , becomes the ground state |0o) of an unmeasured harmonic oscillator with 
frequency uq. In analogy with Eq. (|5.1|), the coherent states in presence of a continuous 
measurement are then defined as 



\p'q') 



(5.24) 
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which, in the position representation, becomes 

-1/4 



Wq') 



exp 



1 - ^(t2 



{q-q') +^p'{<i-(i' 



(5.25) 



with 0"^ and 0"^^ given by Eqs. (p.20|) and ( p.21|) , respectively. The states \p'q') of Eq. 
( p.24| ) have the same properties of the usual coherent states [^. In particular, they form 
an overcomplete basis with completeness relationship 



dp'dq' 
27Th 



\pq'){p'q'\ 



and overlaps 



{p'q \p'q") = exp 



i p+p , „ 



(5.26) 



(5.27) 



where 



Cp'-p"q'-q" 



2h^ 



ll\2 



(5.28) 



Suppose that a quantum system, not necessarily a linear one, is prepared at time t in 
the coherent state \p{t)q{t)). To the leading order in we have 



and 



We also have 



a. 



pq 



(5.29) 



(5.30) 



= {pit)q{tW\p{t)qit)) - {p{t)q{t)\p\p{t)q{t)Y = {<P\p' 

-2 



4 ^ ^ pq 



(Tf, 



(5.31) 



Note that these expressions are independent of ujq 0]. Since a^, a"^, and vanish for h ^ 0, 
in this limit the expectation values p(t) = {p(t)q(t)\p\p(t)q(t)) and g(t) = {pit)q{t)\q\p{t)q{t)) 
can be interpreted as classical phase space coordinates. In a selective measurement, their 
change is given by Eqs. (|5.4|) and ( |5.5D 



dp{t) = - ^p{t) + {p{t)q{t)\dqV{q, t) \p{t)q{t))\ dt + pm^kBTdw^{t) (5.32) 

dw!-{t). (5.33) 



dq{t) = ^dt 



m 



\ SmkBT 
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For ^ — > 0, the expectation value {p{t)q{t)\dgV{q,t)\p{t)q(t)) can be replaced with 
dqV{q(t),t) and we recover the classical Langevin equations ( |2.12| ) and ( p.l3| ). 

In the case of a nonselective measurement, the h ^ limit is properly discussed in 
terms of the Wigner function related to the reduced density matrix operator through the 
transformation 



If f i \ z z 

W{p,q,t) = ^ J dz exp [tVz] {q - i^\Qit)\q + i^) ■ 



(5.34) 



Suppose that at time t the measured system is described by the density matrix obtained by 
averaging the state \p{t)q{t)) over all noise realizations, i.e., all possible values oi pit) and 
q{t) specified by a certain distribution function such that Tr^(t) = 1, 



p(t) = W)q{t)){p{t)q{t)\ = J dp'dq' 6{p' - p{t))6{q' - q{t)) \p'q'){p'q'\. 
The corresponding Wigner function is 



W{p,q,t) = / dp'dq' S{p'-p{t))S{q'-q{t)) Wpy{p,q), 



where 



1 /" ( i \ z z 

Wp>q>{p, q) = ^ J dz exp [rPz] {q - ■^\p'q'){p'q'\q + ^) 



2al 



— exp 
nn 



H 2 



{p-p') - -^{q-q' 



2al 



{q-q 



l\2 



Since 



(5.35) 



(5.36) 



(5.37) 



lini pVpy(P' q) = ^{p- v')5{q - q') 



(5.38) 



in the ^ — limit W{p^ q, t) reduces to the classical probability density 6{p — p{t))6{q — q{t)) 
obtained by averaging the sharp density S{p—p(t))6{q — qit)) over all acceptable phase space 
points p{t), q{t). Finally, the equation of motion of the Wigner function, obtained from Eq. 
( p.36| ) with standard manipulations |^, is 



a,W(p,q,t) 



2n 



m 



n=0 



2i 



+ dp i'jp + m'jkBTdp) + 



W{p,q,t), 



(5.39) 



so that, in the h limit, the change of W{p,q,t) coincides with that prescribed by the 
classical Fokker-Plank Eq. (|2.14| ). 
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VI. MEASUREMENTS ON MACROSCOPIC SYSTEMS 



One of the principal drawbacks of the von Neumann measurement theory is the impos- 
sibihty of predicting a quantum-to-classical transition in the macroscopic limit, unless the 
state \ip{t')) of the system at the beginning of the measurement is one of the coherent states 
( p. 11 ). On the other hand, when the size of the system is sufficiently large, i.e., in the for- 
mal h —>■ limit, we must always recover the result of a classical measurement. It is now 
clearly established that the entanglement of the measured system with the infinitely many 
degrees of freedom of the measurement apparatus can provide superselection rules which 
avoid paradoxical quantum features at macroscopic level 

Concerning the measurement model discussed here, the existence of a superselection rule 
of this kind can be demonstrated in a general way in the case of linear systems. Halliwell 
and Zoupas have shown, in a statistical sense first [jl9[ and with a more direct approach 
but in the free particle case and dissipationless limit later |3^, that the solutions of Eq. 



( p.37| ) converge to a coherent state characterized by time-dependent parameters p{t) and 
q{t) which are the expectation values of p and q in the state itself. Here, we generalize the 
result of by showing that the solutions of Eq. ( |3.37D with potential 

V = Vo{t) + vi{t)q + \mulq^ (6.1) 

in the long time limit become of the form 



= exp 



\p{t)q{t)) (6.2) 



where p{t) = {ipi^]{t)\p\i;i^]{t)) , q{t) = (%](t)|g|V^K](t)) and (fi{t) and \p{t)q{t)) are given by 
Eqs. ( p.23| ) and ( |5.24| ), respectively. 



The Green function corresponding to the nonlinear Eq. (|3.37|) 

G[^]{q, t; q', t') = J d[p]d[q]f^t, exp [j^j^^ds 

+iht,{q - {q) f - ih^{q - (q))^ + TTT^P ' ] (6-3) 



2 

P 1 



functionally depends on the state \ip[s^]{t)) through the expectation values (p) = 

(^) — (V'[C](^)I^I^K](^))- If we suppose, for the moment, that these 
functions and the noise ^ are given, for V of the form ( |6.1|) the Green function ( |6.3|) is that 
of a linear system with Lagrangian 

L{(1, Q, t) = - ^mu'^q'^ - m'jqq + f{t)q + g{t)q + h{t), (6.4) 

where 

m 

and fit), git), and /i(t) are given in terms of (p), (g), ^, Vq, and Vi. By performing the 
Gaussian functional integrals in Eq. (|6.3|), we get 
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Gi^]{q, t] q', t') = n{t, t') exp <j - 



9{t)q - 9{t')q' - + + S^\{q, t] q', t') 



(6.6) 



where n{t,t') does not depend on the spatial variables and 

1 



Sci{q,t; q',t' 



mio 



cos[co'(t — t')] [q + q 



12 



mcuqq 



sm[uj{t -t')][ 

+ qf ds [f{s)-g{s)]sm[Lj{s-t')] + q'f ds [f{s)-g{s)]sm[Lj{t-s)] 
Jt' Jt' 

^ [/(s)-^(s)]sin[a;(t-s)]J^'ciM [/(«)- ()(n)]sin[^(n-t)] I (6.7) 



muj 



is the classical action of a driven harmonic oscillator of mass m, frequency u and external 
force f — g, evaluated with boundary conditions q(t') = q' and q{t) = q [Q. The frequency 
uj is complex with real and imaginary parts given by 



Re(cj) = ± 



h 



h 



2"^ ^ 72 - a;2 + ^(72 - ooiy + {2hK/mf 2ma2 



(6.8) 



Im(ti;) = =F 



\2mal) 



(6.9) 



For {t — t')|Im(a;)| ^ 1, the coefficient of the qq' term in the action ( p.7|) vanishes while the 
coefficient of (g^ + g'^) becomes -^imuj/2. In the long time limit, therefore, the propagator 
( |6.6| ) looses memory of the initial conditions and the solutions of Eq. ( p.37| ) can be written 

as 



V'[e](9,*) = / dq GiQ{q,t;q',t') ij{q',t') 



exp 
exp 



m 



-^(±^ + ^7) q^ + ait)q + Pit) 



1 - #cr2 



h pg „2 



4a2 



q' + a{t)q + I3{t) 



(6.10) 



where we used 



hlm{uj) 

2Re{uj) ' 2Re{uj 



'm2(72-^2) 4 + 



1 2 
m^ya = cr. 



pq- 



(6.11) 



The complex functions a{t) and /3(t) are to be determined. First of all, the normalization 
of the wavefunction ( |6.1(]| ) implies that 



(27raJ)^/2 |2Re[/?(t)] + 2(TjRe[a(t)]^} . 



(6.12) 



Then, we can impose two self-consistency conditions involving the expectation values of p 
and q in the state ( |6.1CI| ) 
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P{t) = {i^amWmit)) = ^Im[a(t)] + 2al^Ma{t)] (6.13) 

and 

qit) = (^K](t)|g|^K](t)) = 2ajRe[«(t)]. (6.14) 

By using the expressions so obtained for Re[a(t)], Im[a(t)], and Re[/5(t)], the wavefunction 
( |6.10| ) can be rewritten as 



1 - fc^^ 



I 



V'[^](g,t) = (2vraJ)-V4exp ^ [q - q{t)f + -p{t)[q - q{t)] - -^{t) \ , (6.15) 



where 



V{t) = -hlm[P{t)] - p{t)q{t) + ^q{tr. (6.16) 

The action if{t) evolves according to an equation obtained by imposing that the change of 
ip[^]{q,t) in a time dt is given by Eq. ( p.37|) . Since the wavefunction (|6.15| ) is of the form 



[2]), the differential dip{t) is given by Eq. ( p.23| ). This completes the convergence proof. 



To the leading order in h, the characteristic time which determines the convergence of 
'0[^](g,t) to the wavefunction ( p.l5| ) is given by 



1 h 



The convergence becomes infinitely fast for h ^ 0. On the base of this result and of the 
properties of the coherent states discussed in the previous Section, we can conclude that 
during a measurement, selective or nonselective, the h —>■ limit does exist at any time 
t > t' even if it does not exist at t = t'. As an example of this behavior, in Appendix B 
we explicitly evaluate the ^ ^ limit in the case of nonselective measurements on a free 
particle cat state. The discontinuity at t = t' is, of course, an artifact of the instantaneous 
correlation assumed through Eq. ( p.3|) between the measured system and the measurement 
apparatus and it would disappear in a more physical approach in which such correlation is 
established in a finite time. 

In the case of nonlinear systems, terms higher than quadratic appear in the potential of 
the Lagrangian ( |6.4| ) so that the convergence proof given for linear systems does not apply. 
However, due to the linearity of the interaction with the infinitely many oscillators of the 
measurement apparatus, the leading h term of this potential, i.e., 

— ihKq^ = —i — (6.18) 

n 



is a quadratic one with complex frequency ^ —Ai^ykBT/h. As a consequence, in the — > 
limit the state of the system becomes of the form (|6.15|) with cr, and apg given by Eqs. (^.29|) 



and ( ^.30|) . The recovering of the classical behavior in the macroscopic limit is, therefore, 
obtained independently of the nature of the measured system. Numerical examples of this 



result can be found in p4|-pq] and experimental evidence has been recently reported in 

After the completion of this paper, we became aware of a preprint by Strunz and Per- 
cival ||3^ in which the authors discuss the semiclassical behavior of open quantum systems 



described by a general Lindblad master equation. 
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APPENDIX A: POSITIVITY AND MARKOVIAN EVOLUTION OF g 

Violations of the positivity of g{t) may arise due to inappropriate approximations of 
the exact influence functional (|3.12| ). Examples of exact solution of g{t) with no positivity 



violations have been given in |^ in the case of a harmonic oscillator. 

In the framework of a theory of measurement processes, the requirement that the meter 
has classical behavior with respect to the observer imposes hfl <C fc^T, if the frequency 
distribution (|2.1CI| ) is assumed. In this case, the fluctuation kernel 



A / \ "^7 f 7 , f \ r / \n 

A( s — u) = — — / au LJ coth — cos c<j(s — u)\ 

Tvh Jo \ 2A;rT / 



has a high temperature expansion whose leading term VksT /h^ is proportional to the dis- 
sipation kernel 

r(t - s) ~ 2m7(5(t - s) (A2) 

which is Markovian for <^t,t being the fastest time scale of the classical motion. When 
these approximations are made in Eq. ( |3.12|) , so that Eq. (|3.15| ) and the corresponding 



master equation (|3.18| ) are obtained, violations of positivity of Q{t) may occur as in the 
example pointed out in However, this happens on a time scale shorter than h/kBT, 
i.e., outside the range of validity r ^ 3> h/ksT of Eq. (|3.18|) []39| . In this range. 



the substantial positivity of Q{t) can be made apparent by selecting appropriate dominant 
terms. 

Equation ( |3.15D shows that the dissipation term is negligible in comparison to the fluc- 



tuation one when 

where qa = — 12 and = Pi + P2- The functions and may assume any value 
according to the functional measure However, close to the dominant classical path 

we have ~ 2mgA/T and the condition ( |A3[ ) can be restated as kBT/h ^ t~^. Therefore, 
in the working range r 3> ^ h/ksT dissipation can be neglected with respect to 
fluctuation and Eq. (|3.15|) becomes 



Z[q,,q,] = I'ds ^^^^ Ms) - q2{s)f . (A4) 
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Correspondingly, the nonselective measurement processes are described by the master equa- 
tion 



rri'ykBT 



(A5) 



This is a Markovian evolution of Lindblad class and therefore (completely) positive |^ . The 
associated selective processes are described in terms of a single state satisfying the stochastic 
Schrodinger equation 



d\^[dt)) = -lH{p,q,t)\ij[^it)) dt 



'^-^^^^^[q-qmi^\imdw,{tl (A6) 



with q{t) = {ilj[(^](t)\q\ip[^](t)) . The general results of are recovered by setting 

At classical level, the Fokker-Plank equation ( |2.14|) and the Langevin equations (|2.12|) 
and ( |2.13| ) are consistent with the fluctuation-dissipation theorem Equations ( |A5| ) 

and ([A^), in which dissipation is neglected, are therefore not appropriate for recovering 
the classical limit. New quantum equations are to be introduced which include dissipation 
and, at the same time, guarantee the positivity of Q{t). As shown in Section III, this is 
accomplished by rewriting equation ( |3.15|) in the equivalent form (|3.21|) and neglecting the 
B terms with respect to the A ones on the base of the working condition h/ksT <^ r. 
Dissipation is still contained in the remaining influence functional which gives rise to the 
master equation ( p.36| ) of Lindblad class and to the corresponding stochastic Schrodinger 
equation ( 3.37| ). These equations provide the correct classical limit as shown in Section V. 

We conclude with some remarks about the possibility pointed out in PT],H2| of obtaining 
a master equation of Lindblad class by taking into account the next to leading term in Eq. 
Jj). In this case, Eq. ( p.l5|) would become 



Z\pi,qi^P2,(l2] = I ds 
Jt' 



- / ds 



du 



UksT 



iQiis) - q2{s)] f (s - u) [qi{u) - q2{u)] + 



(A7) 



The new term can be more easily analyzed after integration by parts. By setting = Q1 — Q2, 
we have 



ds / du qA{s)T(s — u)qA{u) 
Jt' 

= \-fds f du qA{s)T{s - u)qA{u) - \v{^) W{tf + gA(t')'l + r(t - t')9A(t)gA(t') 

I Jt' Jt' 2 L J 

\ ft 1 /■* , , 

+- / du r(t - u) \qA{t)qA{u) - qA{u)qA{t)\ + - / du V{i! - u) \qi\{u)qA{i!) - qA(i!)qA(u)\ . 
2 Jt' 2 Jt' 

For t — f ^ Q^'^, the last three terms can be neglected and the first one approximated with 
a single integral 
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1 /■* 1 1 1 /■* 

- J^^ ds gA(s)gA(s) = -qA{t)qA{t) - ^qAit')qAit') - ^ J^^ ds qA( 



qA{tf + qA{t 



By using the identity 

1 r 
2 

Eq. (|A8|) can be rewritten as |^| 

ds 



l\2 



q^it ) +7; ds gA(s)gA(s) 



(A9) 



(AlO) 



du qAis)T{s - u)qA{u) 
2m7 [qA{t)qA{t) - qA{t')qA{t')] - \ j'^ ds qAi 



-qAitr - - 



j'^^ds gA(s)gA(s)|. (All) 



In Refs. 1 41. 52] the first term of Eq. ( |A1 1|) is neglected by observing that it is much smaller 
than the third one. In this case the reduced density matrix operator git) would undergo a 
transient change 



^(91,92,0) exp 



Mt')-q2it')f\g{q^,q2,0) 



(A12) 



followed by a Lindblad evolution described by an equation which reduces to Eq. ( ^.36[ ) in 
the h ^ and T ^ 00 limits. The validity of these findings is, however, questionable. The 
surface terms neglected in Eq. (|A11|) are of the same order of the integral of ^q\ which is, 
conversely, maintained. 

The existence of a transient change in the evolution of the reduced density matrix de- 
serves further comments. If the system and the measurement apparatus are initially non- 
interacting, a change of g{t) at the switching-on of the interaction is plausible. However, 
as we explained in Section II, this transient cannot be described in the framework of a 
model, such as the bath of harmonic oscillators, in which the system and the measurement 
apparatus are always in interaction. We must limit our considerations to a non-transient 
evolution and, correspondingly, assume that the system and the meter are correlated from 
the beginning. 

The drawbacks of Refs. |41| , |42[| have been recently underlined also in [Q]. Similarly to 
Ref. |2^ and the present work, the authors of assume an initial correlation between the 
system and the environment. 



APPENDIX B: NONSELECTIVE MEASUREMENTS ON A FREE PARTICLE 

CAT STATE 



Let us consider a free quantum particle that, at the beginning of the measurement, is in 
the superposition (cat) state 



mt')) = N Ipigi) + \p2q2) 



(Bl) 



where \piqi), i = 1,2, are two coherent states ( ^.24|) with uq = and the normalization 
factor is 
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The initial state 
Wigner function 



N = 

^2[l + exp(-Cp,_p,g^_,J] 
T|) has no classical counterpart. The h 



(B2) 



limit of the corresponding 



W{p, q, t') = Wp,q, {p, q) + Wp^q^ij), q) + W v,+v,, {p, q) 



x2 cos 



P( N V1-P2 ( qi + q2 



(B3) 



where Wpigi{p,q) is given by Eq. ( |5.37| ), does not exist. The situation changes during the 
measurement. In the case of a nonselective process, the system is described by a Wigner 
function W{p,q,t) which can be exactly evaluated by solving Eq. ( |5.39| ) with the initial 
condition (p3|) [^ . The result is 



W{p, q, t) = N'{ Wp,g,{p, q, t) + Wp,g,{p, q, t) + Wp,+,, ,,+,, (p, q, t)e-^pi-P2«i-.2+^Pi 



x2 cos 



Pi +P2 

2h 



{qi - 92) + T, 



pi-P2qi-q2\^ ) 2 



2m7 



(t\ia- + - ^^+^^ (l - e-^(*-*') 



where 

Wp'q'{p,q,t) 



2ti^ ACUt)Cyy{t) - a,y{ty 

Cxx (t) 



X exp 



Cxy(t) 



m7 ^ 



n 2 



-7(i-t') 



4:Cxx(t)Cyy(t) — C xy (t)'^ 

^yyjt) 

4:Cxx{t)Cyy(t) — Cxy{t)'^ 



q — q ( 1 — e 



-7(t-t') 



(B4) 



(B5) 



4:Cxx{t)Cyy{t) — C xy {t)'^ 



(B6) 



and 



2Cyy{t)Cp, ,{t) — Cxy{t)Cp, ,{t) 



2 ' 



4:Cxx(t)Cyy(t) — C xyit) 



'^Cxxit)Cp,g,{t) — Cxy{i)Cp,q,{t) 

^Cxxif)Cyy{t) — Cxyi^t)"^ 



(B7) 



(B8) 
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are given in terms of the coefficients 



- I - + — 

2m7cr2 I 4 



2 h-f 



Cxyit) — h 



1 




4) 









Cyyit) 



m7 I 2 ^ 



1 - e-^(*-*') 
1 _ 2e-^(*-*') + e-2^(*-*')" 



1-e 

h 



2 r 

g 



1 + ^ 



1 _ e-7(*-0 



+ 



1 ^ 



m7cr^ 



7(t-t')-- + 2e-^(*-*')--e"2^(*-*') 



9 



V4 ■ 



2 



1 ;i7 

16 



-7{t-t') 



7(t - t') 



m70"g 



m7(jg \ 4 



1 - e-^(*-*') 
1 _ e-7(*-*') 



+ 



cr. 



(B9) 



(BIO) 



(Bll) 
(B12) 



(B13) 



The indices x and y stand for the Fourier variables conjugated to p and g, respectively. Each 
term of Eq. (|B4| ) is localized within a phase space region whose size grows with time. When 
this size has became much larger than CTpCTq ~ we can write 



j dp'dq' W{p',q',t) Wp>,ip,q) ~ W{p,q,t) 



(B14) 



and identify the weights 5{p — p{t))5{q — q{t)) of Eq. ( ^.36| ) with W{p,q,t) in agreement 
with |2|. 

Let us now turn the classical limit of Eq. (p^). First, we note that Wp'q'{p,q,t) is the 
time evolution of the Wigner function (|5.371 ) corresponding to the coherent state \p'q')- Its 
classical limit exists and is given by an expression W^'q'{p,(l,t) identical to Eq. (|B5D with 
Cxxit), Cxyit), and Cyy(t) replaced with 



Ct{t) = \iu.Cxx{t) = -mkBT 
h—^o z 



C^lit) = ]im Cyyit) 



7 



1 _ g-27(t-t') 



7 



it-t')-- + 2e"^(*-*') 



-27(t-t') 



(B15) 
(B16) 

(BIT) 



The function Wp}g,ip,q,t) is the phase space probability density obtained by solving the 
Fokker-Plank Eq. ( p.l4| ) with initial condition Wpl^iip, q, t') = 6ip — p')Siq — q'). 

Concerning the interference term in Eq. (p^), we have = Cp>q>, Tp/ = q'/h, 

and $p'g'(t') = —p'/h so that, as previously noted, the h ^ limit of this term does not exist 
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aX t = t' due to the undamped oscillation of the cosine. On the other hand, for t > t' since 
Sp'g'(t) = 0{h~^) while Cpiqi = 0{h~^^'^) we have an exponentially damping term which 
allows to obtain 



1 



Ihn W{p, q,t) = - [W;l^^ (p, g, t) + W;^^^ (p, g, t) 



(B18) 



From a physical point of view, this limit is equivalent to a macroscopic one in which \pi — 
V'AI'^p and/or — Oq become infinitely large so that Cpj.pggi-ga diverges. In particular, 
this is obtained by taking the mass m of the particle infinitely large. 

Finally, we note that due to the condition -C k^T the Wigner function Vrp/g/(p, g, t) 
approaches the classical phase space probability density W'^^iip.q^^) on a time scale 
{h/'^kBTY^'^ <^ On the other hand, the functions Sp/g'(t), Tp/q/(t), and $p'g'(t) 

vanish for t —>■ oo. The long time limit of Eq. ([B4|), therefore, is 



xe '^pi-P2M-g2 2cos 



P1+P2 
2h 



(B19) 



and never coincides with the classical limit ( [B18|) . 
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